ICTP-SISSA 



Coherent oscillations between two weakly coupled Bose-Einstein 
Condensates: Josephson effects, 7r-oscillations, and macroscopic 

quantum self trapping 

S. Raghavan 1 , A. Smerzi 2 , S. Fantoni 1 ' 2 , and S. R. Shenoy 1 
1) International Centre for Theoretical Physics, 1-34100, Trieste, Italy 
2) Istituto Nazionale de Fisica della Materia and International School for Advanced Studies, via 

Beirut 2/4, 1-34014, Trieste, Italy 
(February 1, 2008) 

Abstract 



We discuss the coherent atomic oscillations between two weakly coupled 
Bose-Einstein condensates. The weak link is provided by a laser barrier in a 
(possibly asymmetric) double-well trap or by Raman coupling between two 
condensates in different hyperfine levels. The Boson Josephson Junction 
(BJJ) dynamics is described by the two-mode non-linear Gross-Pitaevskii 
equation, that is solved analytically in terms of elliptic functions. The BJJ, 
being a neutral, isolated system, allows the investigations of new dynamical 
regimes for the phase difference across the junction and for the population 
imbalance, that are not accessible with Superconductor Josephson Junctions 
(SJJ). These include oscillations with either, or both of the following prop- 
erties: 1) the time-averaged value of the phase is equal to ir (it — phase 
oscillations); 2) the average population imbalance is nonzero, in states with 
"macroscopic quantum self-trapping" (MQST). The (non-sinusoidal) gener- 
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alization of the SJJ 'ac' and 'plasma' oscillations and the Shapiro resonance 
can also be observed. We predict the collapse of experimental data (corre- 
sponding to different trap geometries and total number of condensate atoms) 
onto a single universal curve, for the inverse period of oscillations. 

Analogies with Josephson oscillations between two weakly coupled reser- 
voirs of 3 He-B and the internal Josephson effect in 3 He- A are also discussed. 
PACS: 03.75.Fi,74.50.+r,05.30.Jp,32.80.Pj 
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I. INTRODUCTION 



Bose-Einstein condensation, predicted more than 70 years ago 0, was detected in 1995, 
in a weakly interacting gas of alkali atoms held in magnetic traps 0. Following the first 
observations, there have been important experimental developments. A superposition of 
condensate atoms in different hyperfine levels has been created; non- destructive, in situ, 
detection probes have tracked the dynamical evolution of a single condensate ||. Further, the 
evolution of the relative phase of two condensates has been measured through interferometry 
techniques 0. More recently, experiments that tune the scattering length by several orders 
of magnitude have opened the definite possibility of creating in the laboratory an ideal 
condensate of non-interacting atoms. 

The precise manipulation of this form of matter is of considerable theoretical interest: 
besides the study of fundamental aspects of superfluidity from "first principles" , it is possible 
to address "foundational" problems of quantum mechanics ||. In fact, the order parame- 
ter can be identified with the one-body macroscopic condensate wave function. This obeys 
a nonlinear Schrodinger equation, known in the literature as the Gross-Pitaevskii equation 
(GPE) JJ. The GPE has been successfully applied to study kinetic properties of the conden- 
sate, like collective mode frequencies of trapped Bose-Einstein condensates (BEC) [Kj and 
the relaxation times of monopolar oscillations |TT| . Chaotic behavior in dynamical quantum 
observables [|TT1 , |l"2f , and the metastability of quantized vortices has been predicted . 

The existence of spatial quantum coherence was demonstrated by cutting a single trapped 
condensate by a far off-resonant laser sheet. Upon switch-off of the confining trap and the 



laser, the two condensates overlapped, producing interference fringes |[L4|| , in analogy with 
the a double-slit experiment with single electrons. 

However, the superfluid nature of BEC can be fully tested only trough the existence of 
superflows. Current experimental efforts are being focused on the creation of a Josephson 



junction between two condensate bulks fll4yi5| . In this context, the Josephson junction 
problem has been studied theoretically in the limit of non-interacting atoms ||16||, for small 
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amplitude Josephson oscillations [|T7||l"8|l , including finite temperature (damping) effects |L8 . 
Furthermore, in connection with the quantum measurement problem, decoherence effects 
and quantum corrections to the semiclassical mean-field dynamics have also been 



studied. Self trapping dynamics in the limit of small number of condensate atoms has been 
considered |]TU| in the "quantum" and in the "semiclassical" (mean-field) approximation. 



We have elsewhere, [21]], pointed out that even though the Boson Josephson Junction (BJJ) 
is a neutral-atom system, it can still display the (non-sinusoidal generalization of) typical 
'dc', 'ac' and Shapiro effects occuring in charged Cooper-pair superconducting junctions. 
Moreover, novel dynamical regimes like macroscopic quantum self trapping (for arbitrarily 
large condensates) and 7r-phase oscillations (where the average value of the phase across the 
junction is equal to it), have been predicted. In the present paper we present a comprehensive 
analysis of the effects described in pl| , including a discussion of the BJJ equations and 



their analytic solution, and limits of the approximations underlying the BJJ model, and a 
comparison with other superconducting and superfluid Josephson junctions. 

The description of the GPE dynamics for a Bose condensate in a double-well trap re- 
duces, under certain conditions, to a nonlinear, two-mode equation for the time- dependent 



amplitudes ^1,2 (*) = \J N li2 (t)e i()l - 2{t \ where N h2 (t) and 6 lj2 {t) are the number of atoms and 
the phases of the condensate in the trap 1, 2 respectively. These amplitudes are coupled 
by a tunneling matrix element between the two traps, with the spatial dependence of the 
GPE wave function integrated out into constant parameters. The resulting BJJ tunneling 
equations resemble the (nonlinear generalizations of) Superconductor Josephson Junction 
(SJJ) equations, with the variables being the relative phase and the fractional population 
imbalance. 

However, there are important physical differences between the isolated double-well BJJ 
and the SJJ with an external circuit. The SJJ is generally discussed in terms of a rigid 
pendulum analogy in the resistively shunted junction (RSJ) model, while the BJJ in a double- 
well trap can only be completely understood in terms of a non-rigid pendulum analogy, 
with a length dependent on the angular momentum. In SJJ the Cooper-pair population 



imbalance is identically zero (considering two equal-volume superconducting grains) due to 
the presence of the external circuit |22| , and the dynamical variable is the voltage ~ </> across 
a quasiparticle resistive shunt. In the BJJ, the non-rigid pendulum dynamics are associated 
with superfluid density oscillations, of an isolated system. An isolated (without external 
circuit) superconducting junction allows coherent Cooper-pair oscillations, but only in the 
small amplitude (plasma) limit p2[-|24j . 

A closer analog of the BJJ is provided by the internal Josephson effect in 3 He-A, where 
the (rigid) pendulum oscillations describe the rate of change of up and down spin-pair 



populations, induced by an external variable magnetic field p5|-p7||. The u ir oscillations" 
between two weakly coupled reservoirs of 3 He-B [28] could be related to the analogous 
oscillations occuring in BJJ. 

The experimental detection of predicted effects in BJJ could be achieved through tempo- 
ral modulations of phase-contrast fringes ||14| , interferometric techniques ||, or other probes 
of atomic populations , using ~ millisecond temporal oscillations of the integrated sig- 
nal Ni — N 2 . The direct detection of the currents instead of densities, perhaps by Doppler 
interferometry, would be worth exploring. 

The plan of the paper is as follows: In Section |TJ, we obtain the BJJ tunneling equations, 
that, in Section |Tl|, are compared with the Josephson equation for other superconductor 
and superfluid systems. In Section fV| we solve the BJJ equations discussing the various 
dynamical regimes. In Appendix |A| we outline the derivation of the two-mode BJJ from 
GPE and discuss the limit of the approximations. The BJJ equations are solved analytically 
in terms of elliptical functions in Appendix [FJ. In the remainder of the paper Section [V], we 
discuss the asymmetric trap case, clarifying the analogies with the ac and Shapiro effect. 
We summarize our results in Section [VI. 
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II. BOSON JOSEPHSON JUNCTION: THE NON-LINEAR TWO-MODE 

APPROXIMATION 



The wavefunction ^(r, t) for an interacting BEC in a trap potential V trap (r, t) at T = 
satisfies the GPE: 

ih ^A = _A. V 2 lIr(rt) + [V trap (r)+g \^(r,t)\ 2 }^(r,t), (2.1) 

with go = 4:7ih 2 a/m, m the atomic mass, and a the s-wave scattering length of the atoms 
In the following we will consider a double- well trap produced, for example, by a far 
off-resonant laser barrier that cuts a single trapped condensate into two (possibly asymmet- 
ric) parts ||14|| . However, the results could also apply to the oscillations of the condensate 



population difference between two hyperfine levels fllSf . 

Since we are interested in the dynamical oscillations of the two weakly linked BEC, we 
write a (time-dependent) variational wave-function as: 

tt(r,t) =ijj x {t)<& x {r)+i) 2 {t)<$> 2 {r) (2.2) 

with ipi t2 (t) = \J ' Ni^e 101 ' 2 ^ and a constant total number of atoms N\ + N 2 = |^i| 2 + \ip 2 \ 2 = 
Nt- The amplitudes for general occupations iV li2 (t), and phases 12 (t) obey the nonlinear 
two- mode dynamical equations [|T^-^T],^T],^[ . 

ih^ = {E° x + UxN^fa - /C^ 2 (2.3a) 
ih^- = (E° 2 + U 2 N 2 )i, 2 - Kfa, (2.3b) 

where damping and finite temperature effects are ignored. Here E® 2 are the zero-point en- 
ergies in each well, U\ i2 N\ i2 are proportional to the atomic self-interaction energies, and /C 
describes the amplitude of the tunneling between condensates, see Fig. |I[ The constant pa- 
rameters E± 2 , {Ti.,2, can be written in terms of $>i t2 (r) wave-function overlaps. The &i, 2 (r), 
roughly describing the condensate in each trap, can be expressed in terms of stationary 
symmetric and antisymmetric eigenstates of GPE, (see Appendix |A]). 
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The fractional population imbalance 

z(t) = (N^t) - N 2 (t))/N T = (|^| 2 - \H 2 )/N T (2.4) 

and relative phase 

( j ) (t) = 9 2 (t)-9 1 (t) (2.5) 

obey 

z(t) = -y/l - z 2 (t) sin(0(t)), (2.6a) 

0(t) = AE + Az(t) H t= cos(0(t)), (2.6b) 

'1 - z 2 (t) 



where we have rescaled to a dimensionless time, t2K/% — > t, and 

ee (£° - ^)/2/C + (C/i - f/ 2 )/2/C, (2.7a) 
A ee UN T /2K ; U=(U 1 + U 2 )/2 (2.7b) 

The dimensionless parameters A and AE determine the dynamic regimes of the BEC atomic 
tunneling. The total, conserved energy is: 



Az 



2 



H = — — h AEz — Vl — 2 2 cos 0. (2.8) 



(2.9) 



suggesting that the equations of motion (|2.6|) can be written in the Hamiltonian form 

dH . d# 



dcj) dz 
with z and 0, the canonically conjugate variables. 

For well defined mean values in relativ population and phase, fluctuations must be small. 

III. THE JOSEPHSON EFFECT IN OTHER SUPERFLUID AND 
SUPERCONDUCTING SYSTEMS 



The superconducting Josephson Junction (SJJ). We now consider the SJJ dy- 
namic equations p2|- |2^j33[| , for comparison with the BJJ tunneling equations of Eq. fl2.6|) . 
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The SJJ has an external closed circuit that typically includes a current drive I ex t; the mea- 
surable developed voltage across the junction V is proportional to the rate of change of 
phase: 

dV V 
hxt = Cj— + Ijsmcj) + — (3.1a) 

(3.1b) 

where Cj(Ij) is the junction capacitance (critical current) and R is the effective resistance 



offered by the quasiparticle junction and the circuit shunt resistor. The yl — z 2 factors 
of Eq. ( |2.6| ) are missing here, since the external circuit suppresses charge imbalances, i.e. 
z(t) = The junction charging energy Eq ~ Cj 1 ; superconductor-grain charging 

energies Ecg{~ inverse grain sizes) that are the analogs of interatom interactions U of the 



BJJ, are relevant only in mesoscopic systems. Two such small isolated grains [34| can be 
considered a closer superconducting analog of BJJ. Even in that CclSC, clS Nt is still large, 
the voltages that appear are 2eV ~ 2A gp , the quasiparticle gap, implying that \z\ ~ 10~ 9 . 
Mechanical analogs have been useful in visualizing the SJJ. Equation (|3.1| ) can be written 

as 

4> + ^/RCj + ouj sin = (I ext /Ij)u 2 j (3.2) 



in unsealed units, with uj = y/E c Ej/h, the Josephson plasma frequency. This can be 
regarded as the equation for a 'particle' of 'mass' ~ uj 2 and 'position' <p moving on a 
tilted, rigid 'washboard' potential — cos0— (I ext /Ij)<f), with friction coefficient ~ 1/RCj. 
Alternatively, Eq. ( |3.1|) describes [|33| a rigid pendulum of tilt-angle 0; moment of inertia 



~ ujj 2 ] angular momentum V oc the angular velocity; damping rate (RCj) 1 and external 
torque ~ I ext . 

The Josephson effects in SJJ follow at once from physical considerations: 
Plasma oscillations: For I ext = 0, the rigid pendulum can have small, harmonic oscillations 
in angle around the vertical. Linearizing Eq. ( ^.1[ ) produces sinusoidal volt age /current 
'plasma' oscillations of angular frequency (in unsealed units) 
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u w lo p = 2n/r p = ^E c Ej/k, (3.3) 
independent of initial conditions 0(0), 0(0). 

ac effect: In the pendulum analogy, the external drive balanced by the damping enforces 
steady rotatory motion for I ext /Ij > 1. The phase increases linearly with time, <j>(t) ~ 
2eVt/h, where V = I ex tR is the dc voltage developed, and the current oscillation has angular 
frequency 

2tt 2eV , . 

u = uj ac = — = — — , (3.4) 

Tac n 



independent of 0(0), 0(0). 

Shapiro resonance effect If a small ac component is added to an applied dc voltage, AE — > 
AE(1 + 5 cosu t)] 5q <C 1, then at resonance uj q = u ac , there is a dc tunneling current with 
a nonzero time average (z(t)) ~ 5o(sm.(u ac t + 0(0)) sinc^ot) 7^ 0. This Shapiro resonance 
repeats at higher harmonics u ac = 27r/r ac = ruvo,n = 1,2,..., with characteristic Bessel 
function coefficients J n (n5 ) |22]j23[ . 



Can the BJJ show the full range of SJJ effects? Not at first sight, since the double-well 
BEC is a neutral-atom system. However, the ability to tailor traps and the condensate self- 
interaction compensates for electrical neutrality Asymmetric positioning of the laser 
barrier could produce a zero-point energy difference AE, analogous to an applied voltage, 
since the effective potential seen by the atoms on the smaller-volume side will have a larger 
curvature. The interwell difference between the (bulk) nonlinear atomic self-interaction, 
~ UNtz plays the role of a junction capacitance energy in the dynamics. 

In the SJJ, Ej as well as iV^ are fixed [P^ , ^| . For the BJJ, the laser-sheet intensity, 
and hence the coupling /C, can be varied. Initial states, iVi(0) 7^ A^(0), i.e., z(0) 7^ 0, can 
be prepared, and the laser barrier then lowered to permit tunneling. 

The internal Josephson effect in 3 He-A. A closer analog of the BJJ Eqs. (|2.6| ) is 



provided by the longitudinal magnetic resonance in He-A [p5|1 , that is generally understood 



as internal Josephson oscillations between two interpenetrating populations of superfluid 
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up-down spin-pairs | 26|| . The weak coupling is provided by the dipole interaction between 



pairs of up and down spins. The spin dynamics is governed by p7[ : 



z(t) = -sinO(t)), (3.5a) 
<p(t) = AE + Az(t) (3.5b) 

where z(t) is the fractional population imbalance between up and down spin Cooper 
pairs, A oc (xQd)' 1 with x-,9d the susceptibility and the dipole coupling respectively, and 
AE oc (B/x9d) with B the external applied static magnetic field. In |25| experiments 



have confirmed Eqs. ( |3.5|) , showing the transient between the small amplitude and ringing 
oscillations of the pendulum equations Eqs. ( |3.5| ). 

Josephson Oscillations between two weakly linked reservoirs of 3 He-B. Quite re- 
cently the first direct experimental observation of Josephson oscillations between two weakly 
linked superfluid systems has been reported p5| , |36f . The weak link was provided by ~ 4000 
small holes in the rigid partition separating two 3 He-B superfluid reservoirs, with the hole 
diameter being comparable to the coherence length. A soft membrane created a pressure dif- 
ference across the weak link, inducing Josephson mass current oscillations. These oscillations 
obey: 

I(t) = / c sin(0(*)), (3.6a) 
0(t) = -^AP (3.6b) 

hp 

with 2m3 the mass of a 3 He Cooper pair, p the liquid density and AP the pressure difference 
across the weak link being proportional to the elastic constant of the membrane. Small and 
large amplitude oscillations have been observed, as well as the driven running solutions of 
the phase — oo < <fi < oo [36|| , corresponding to a self-maintained population across the weak 
link. 

By driving the soft membrane in resonance with the natural Josephson frequency, a 
new metastable dynamical regime was observed, with the time-averaged value of the phase- 
difference across the junction equal to 7T. These metastable '7r-oscillations' have amplitudes 
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and frequencies smaller than the 'stable' Josephson oscillations, into which they decay with 
a life time that increases with decreasing temperature Analogous 7r-oscillations with 
similar properties are described by BJJ (see Section |TV| ). In a different context, 7r-junctions 
have been created with high-T c superconductors, that reflect the symmetry of the d-wave 
pairing state 

IV. THE SYMMETRIC TRAP CASE, AE = 

Stationary Solutions. For a symmetric BJJ, i.e., AE = 0, the equations of motion 
Eq. |2]6| are 

z(t) = -y/l-z*{t)sm(<t>(t)), (4.1a) 

kt) = Az(t) + ; = cos(0(t)), (4.1b) 
y/l-z*(t) 

with the conserved energy: 

Az(0) 2 i 

H = H(z(0),<P(0)) = " V 1 " ^(0) 2 cos 0(0). (4.2) 

The ground state solution of the symmetric BJJ, Eq. ( |4.1|) , is a symmetric eigenfunction 
of the GPE with energy E + = — 1 and: 

4> s = 2nn (4.3a) 
z s = (4.3b) 

The next stationary state at higher energy EL = 1, is an antisymmetric eigenfuncion with: 

<p s = (2n + 1)tt (4.4a) 
z s = (4.4b) 

For non-interacting atoms in a symmetric double-well potential, the eigenstates of the 
Schrodinger equation are always symmetric or antisymmetric, with z s = 0. However, be- 
cause of the nonlinear interatomic interaction, there is a class of degenerate GPE eigenstates 
that break the z — symmetry: 
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(p s = (2n + 1)tt (4.5a) 
*. = ±\/l-^ (4.5b) 

provided that |A| > 1. The energy for this state is = | fA + -jm. 

These z-symmetry breaking states are an artifact of the semiclassical limit in which the 
GPE has been derived. In a full quantum two mode approximation the eigenstates are 
always symmetric in the population imbalance: as we will discuss later, such states have a 
large lifetime that scales exponentially with the total number of atoms. 

Rabi Oscillations. For noninteracting atoms (A = 0) Eqs. ( |4.1|) describe sinusoidal 
Rabi oscillations between the two traps with frequency lor = |/C. These oscillations are 
equivalent to single atom dynamics, rather than a Josephson-effect arising from the inter- 
acting superfluid condensate. The possibility to tune the scattering length to values very 
close to zero 0, open avenues for their experimental observation. 

Zero-phase modes. These modes describe the intra-well atomic tunneling dynamics 
with a zero time-average value of the phase across the junction, {<f>(t)) = 0, and (z) = 0. To 
this dynamical class belong small and large amplitude condensate oscillations. 

Small amplitude oscillations. The small amplitude, or "plasma" (in analogy with S J J), 
effect follows at once from the pendulum analogy. From Eq. (|4.1|) , the BJJ is like a non-rigid 
pendulum of length 

(x 2 + y 2 ) 1/2 = (4.6) 

decreasing with angular momentum z, and with moment of inertia A -1 . Linearizing 
Eq. (|4.1|) , we obtain sinusoidal oscillations with inverse periods (in unsealed units) 



r L 1 = J2UN T K, + (2/C) 2 /27rfr, (4.7) 



independent of the initial conditions z(O),0(O). The comparison between Eq. ( |4.7| ) and 
Eq. ( |3.3| ) indicates that 2N T fC (~ N T ) is the analog of the Josephson coupling energy Ej, 
while U (~ N^^ 5 , with d the dimensionality of the system) is the analog of the capacitive 
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energy Ec- Since the coupling energy, fixed by the laser profile, is K, ~ A, the tunnel junction 
area whereas the bulk interaction, UN?, is independent of A, the oscillation rate goes as 
t^ 1 ~ A x l 2 . (The plasma frequency for SJJ, r^ 1 ~ \JE c Ej by contrast, is independent of 
A, as Ej ~ A,E C ~ A" 1 ). 

The Josephson-like length £ Aj' = \Jh 2 /2m/C, that governs spatial variation along the 
junction, should be much greater than ~ v^4 to justify neglect of spatial variations of z, 0, 
i.e., to obtain a 'flat plasmon' spectrum. For /C = O.lnK, one finds 'Aj' ~ 10/im. We will 
not, however, consider such spatial variations here. The frequency of the small amplitude 
oscillations in BJJ are of the order of 10 — 100 Hz for typical trap parameters, and should 
be compared with the frequencies of plasma SJJ that are of the order of GHz. 

Large amplitude oscillations. In Fig. [| we display this regime of anharmonic oscillations, 
plotting z = Nl ^ 2 as a function of time, with the initial value of the phase difference 
0(0) = and A = 10, and for increasing values of the inital population imbalance z(0). 
Specifically, z(0) takes on the values 0.1, 0.5, 0.59, 0.6, and 0.65 from a) through e) respec- 
tively. Increasing z(0) for fixed A, (or increasing A for fixed z(0)) adds higher harmonics to 
the sinusoidal oscillations, corresponding to large amplitude oscillations of the (non-rigid) 
pendulum. This is shown in Fig. |2| (b,c). The period of such oscillations increases with z(0), 
then decreases, undergoing a critical slowing down (Fig. 0(d), dashed line), with a logarith- 
mic divergence. The singularity in the period corresponds to the pendulum in a vertically 
upright position, i.e., reaching the fixed point of Eq. (|4.4b]) . 

Running-phase modes: Macroscopic Quantum Self- Trapping. In addition to 
anharmonic and critically slow oscillations, other striking effects occurin BJJ. For instance, 
for a fixed value of the initial population imbalance, if the self-interaction parameter A 
exceeds a critical value A c , the populations become macroscopically self-trapped with (z) ^ 
0. There are different ways in which this state can be achieved and all of them correspond 
to the condition (which we shall term the MQST condition) that 

H = H(z(0), 0(0)) = ^(0) 2 - ^l-^(O) 2 cos(0(O)) > 1 (4.8) 
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In a series of experiments in which 0(0) and z(0) are kept constant but A is varied (by 
changing the geometry or the total number of condensate atoms, for example), the critical 
parameter for MQST is 



l + v /l-.(O)2cos(0(O)) 

z (oy/2 • 1 yj 

On the other hand, changing the initial value of the population imbalance z(0) with a fixed 
trap geometry and total number of condensate atoms (and initial value 0(0)), A remains 
constant and Eq. (|4.8|) defines a critical population imbalance z c . As we shall see in this and 
the next section, for 0(0) = 0, if z(0) > z c , MQST sets in, but for 0(0) = n, z(0) < z c marks 
the region of MQST. More generally, if |0(O)| < 7r/2, MQST occurs for z(0) > z c while for 
other values of 0(0), it occurs for z(0) < z c . 

In this section, we will discuss the type of MQST in which the phase difference of the 
order parameter across the BJJ runs without bound; other types are discussed later. The 
phenomenon can be understood through the pendulum analogy. If the population imbalances 
are prepared such that the initial 'angular kinetic energy' of the pendulum, z 2 (0), exceeds the 
potential energy barrier height of the vertically displaced = tt 'pendulum orientation', there 
will occur a steady self-sustained 'pendulum rotation', with nonzero angular momentum (z), 
and a closed-loop trajectory around the pendulum support. For H < 1 the population 
imbalance oscillates about a zero value. For H > 1 the time-averaged 'angular momentum' 
is nonzero, (z(t)) ^ 0, with oscillations around this nonzero value (Fig. 0). MQST is a 
nonlinear effect arising from the self-interaction ~ UN^z 2 of the atoms. It is dependent on 
the trap parameters, total number of atoms, and initial conditions, and is self-maintained 
in a closed conservative system without external drives. Although the SJJ 'ac' effect in the 
RCSJ model involves a running-phase, it is clearly physically different from MQST, as it is 
a driven steady-state independent of initial conditions. Moreover, in SJJ the Cooper pair 
population imbalance are locked to zero by the external circuit. MQST differs from single- 
electron Coulomb blockade effect that involves a single electron. It also differs from the 



self-trapping of polarons |32[| that arise from single electrons interacting with a polarizable 
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lattice: arising, instead, from self-interaction of a macro scopically large, number of coherent 
atoms. 

7r-phase modes These modes describe the tunneling dynamics in which the time- 
averaged value of the phase across the junction is ((f)) = n. 

The modes arise once more from the non-rigidity (momentum dependent length) of the 
pendulum and are not observable with SJJ. Thei include small amplitude, large amplitude, 
and macroscopic self-trapped oscillations. The last has nonzero average population imbalanc, 
while < z >= for the others. We summarize this behavior in the temporal evolution of 
z(t) in Fig. | for z(0) = 0.6,0(0) = ir. A takes the values 0.1, 1.1, 1.111, 1.2, 1.25, and 1.3 
in Figs. |3|(a-f) respectively. 

Small amplitude oscillations. For small z, Eqs. ( ^4.1| ) can be linearized around the fixed point 
( |4.4b|) yielding harmonic oscillations for A < 1, with a period (in unsealed units): 



t-- 1 = ^/(2/C) 2 - 2UN T JC/2irh. (4.10) 
It is worth noticing that the ratio of the frequency of the small amplitude zero- and 7r- mode 



phase oscillations is ^ = y y+x < 1 (similar to the 3 He-B 7r-oscillations of Section |TJ). 
Linearizing Eqs. ( flip in z only, the BJJ Eq. ( f4.1b| ) reduces to the very simple form: 

4>=-[A sin(0) + - sin(20)] + 0(z 2 ) (4.11) 

This suggests a mechanical analogy in which a particle of spatial coordinate <fi moves in 
the potential: 

V{4>) = -A cos(0) - - cos(20) + 0{z 2 ) (4.12) 

In Fig. |] we see that V(4>) has a small valley around <fi = tt where the particle can oscillate. 
The depth of this valley decreases as A — > 1. The valley persists, in the full potential for 
V((ft), retaining all the higher order terms in z. 

Large amplitude oscillations For -zr-phase oscillations, the momentum dependent length 
allows the pendulum bob to make inverted anharmonic oscillations with (z) = around the 
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(top of the) vertical axis. For large amplitude z(t) oscillations, A can exceed unity, as shown 
in Fig. |(b). 

oscillations with macroscopic quantum self-trapping Here teh non-rigidity allows the 
penulum bob to make a closed < z >^ rotation loop around the top of the vertical 
axis. There are two kinds of such 7r-phase modes with MQST: those where the time average 
(z) < \z s \ 7^ 0, and those where (z) > \z s \ ^ 0, with z s being the stationary ^-symmetry 
breaking value of the GPE. These two kinds of MQST are shown in the time evolution 
of z(t) in Figs. ||](d-f). In Fig. 0(d), the system is in the first type of trapped state. A 
change over occurs at the stationary state (Fig. ^(e), dashed line). Once A exceeds this value 
A s = 1/ yj\ — z(0) 2 (cf. Eq. |4.5|) , the system goes into the second type of 7r-phase trapped 
state (Fig. |(f)). 

In order to see these different kinds of (running- and 7r-phase) MQST modes more trans- 
parently, one can use the energy H = H of Eq. ( |4.2| ) to write the system of equations ( |4.1|) 
in terms of an equation of motion of a classical particle whose coordinate is z, moving in a 
potential W(z) with total energy Wo, 

zitf + W(z) = W , (4.13a) 

where 

W{z) = z 2 (l- AH + ^-z 2 ^J ; Wo = W(z(0)) + z(0) 2 (4.13b) 

Figure |5| displays the potential W(z) against z (Figs. ||(a) and (c)) and the corresponding 
evolution of 0(i) (Figs. |5| (b) and (d)) to display the various dynamical regimes. In Fig. ||(a) 
and (b) 0(0) = 0,A = 10 and in Fig. § (c) and (d) 0(0) = vr, A = 2.5. The horizontal 
lines indicate the energy value Wq. For fixed value of A, and 0(0) = 0, increasing the value 
of 2i(0) changes W(z) from a parabolic to a double-well. The motion of the particles lies 
within the classical turning points in wich the total energy equals the potential energy. For 
z(0) = 0.1, in Fig. |5](a), the potential is parabolic, and the (small amplitude) oscillations are 
sinusoidal. For z(0) = 0.6 the trajectory of z(t) becomes markedly nonsinusoidal, given the 
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double-well structure of W(z). For z(0) > 0.6 the total energy is smaller than the potential 
barrier, forcing the particle to become localized in one of the two wells. The symmetry of 
the classical motion is broken. This corresponds to an MQST state. Figure |5](b) displays the 
corresponding phase = arccos[(A,2 2 /2 — H )/\/l — z 2 } versus z. For untrapped oscillations, 
the 0, z trajectory is a closed curve, with a time-averaged value of 0(t) = 0. In the running- 
mode MQST regime — oo < 0(t) < oo. for the corresponding evolution. 

Let us now focus our attention on Figs. H(c,d). For A = 2.5,0(0) = ir, the ^-potential 
always has a double- well structure and the system is self-trapped for all values of z(0). 
For small values of z(0), the phase 0(t) is unbounded and the system exhibits running- 
phase MQST. However, above a certain value of z(0) = 2z s = 2yl — 1/A 2 , the phase 0(t) 
becomes localized around tt and remains bounded for all larger values of z(0). In Fig. |5](c,d), 
z(0) = 0.7 and 0.98 mark the two different kinds of 7r-phase MQST since they are on either 
side of the stationary state value of z s = ^1 — 1/A 2 . This point will become more clear in 
the phase-plane portrait of Fig. |7]. 

Discussion of Results. A clear observational feature of the behavior of the system 
is the time-period of oscillations. To this end, we plot in Fig. [], the inverse period 1/r 
as a function of the ratio between the initial population imbalance z(0) and the critical 
population imbalance z c . Figure |6|(a) shows the case for 0(0) = and A = 10 (z c = 0.6) 
(dashed line) and A = 100(z c = 0.2) (solid line). The initial parts of the graph for z(0) <C z c 
mark sinusoidal small amplitude ('plasma') oscillations (Fig. 0(a)). On increasing z(0), the 
oscillations become highly anharmonic, with the inverse period that first increases, then 
decreases, displaying a critical slowing down. The logarithmic divergence of the time-period 
at z(0) = z c is marked by the hyperbolic secant evolution of z(t) (Fig. 0(d)). In the inset, 
we show the average value, (z), as a function of z(0)/z c . MQST is signaled by the sharp 
('phase-transition'-like) rise of (z) from zero to a non-zero value. For 0(0) = it, Fig. |^(b), 
something different happens. MQST occurs for values of the initial imbalance z(0) less than 
z c . At z(0) = z c the time-period diverges, and for larger values of z(0), MQST disappears. 

The dynamical behavior of the BJJ system can be summarized quite conveniently in 
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terms of a phase portrait of the two dynamical variables z and 0, as shown in Fig. |7|. The 
trajectories are calculated for different values of A/A c with z(0) kept constant at 0.6. The 
light solid lines mark the evolution for the evolution where the phase oscillates around 
and (z) = 0. The running mode MQST is shown by the trajectories with small dots for 
A/A c = 1, 1.5 with the initial condition being 0(0) = 0. Note that for a rigid pendulum 
(without the y/1 — z 2 term in the Hamiltonian in Eq. (|2.8|) ), one would obtain only the 
curves described thus far. However, for the BJJ, due to the momentum dependent potential 
in Eq. (|2.8|) , there is considerable richness as exhibited by the dark solid lines, dashed lines, 
and lines with large dots. All these curves correspond to 0(0) = tt. Note, for instance, that 
as A/A c increases and approaches unity, the area enclosed by the trajectory shrinks and is 
pinched at A = A c marking the onset of 7r-phase MQST with < z >< \z s \ (dashed line). 
Upon further increase, the area collapses to a point at the z = z s stationary ^-symmetry 
breaking state. Further increase of A/A c induces a reflection of the trajectory about the fixed 
point, and 7r-phase MQST with < z » \z s \ (lines with large dots). Finally, the trajectories 
join the running mode MQST for A/A c = 2.7 (lines with small dots). 

We now outline a possible procedure for a) bridging experimental data with our the- 
oretical model, and b) collapsing data from different experiments onto a single universal 
curve. We note, at the very outset, that other procedures could be experimentally more 
accessible, particularly since novel methods of tailoring traps || and the possibility of tun- 
ing the scattering length of atoms JF] have become current. The calculation of the values 
of A and K, from the experimental data (for a given trap geometry and total number of 
condensate atoms) is straightforward. The onset of MQST is provided by Eq. ( f4.9| ), which 
immediately gives the value of A c from the (experimentally imposed) initial conditions z(0) 
and 0(0). Moreover, in the small amplitude limit the inverse period of the oscillations, given 
by Eq.(|4.7|), provides the value of K from the previously calculated A. 

Different experiments done by varying arbitrarily the trap geometry and the number of 
condensate atoms give a set of parameters A, /C. The data collapse onto a single universal 
curve of tt/(JCCA 2 t) versus k 2 (A) of Eq. (p7| ), as shown in Fig. |8|. 
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The parameters UNt, and E° can be estimated to be ~ lOOnK and ~ lOnK respectively 
for N T = 10 4 if we take the trap-frequency u) trap to be ~ 100Hz. A = UN T /2K can be 
varied widely by changing Nt or the barrier height ~ K, that depends exponentially on 
the laser-sheet thickness. Typical frequencies are then 1/tl ~ 100Hz. With collective 
mode excitation energies A col [ ~ E°, and quasiparticle gaps A qp ~ y/UNxE , for UNtz < 
A qp ^ co u intra-well excitations are not induced. At nonzero temperatures, BEC depletion and 
thermal fluctuations will renormalize the parameters in Eq. (|2.7f ), and will damp Jl8| the 
coherent oscillations. The effects of damping on the oscillation behaviour requires a separete 
treatment and will be considered elsewhere. 



V. THE ASYMMETRIC TRAP CASE, AE / 

Exact solutions and temporal behavior. Let us now consider the case where the 
traps are asymmetric, i.e., AE ^ 0, as in Fig. [1], with Hamiltonian 

Az 2 , 

H = — + AEz - Vl- z 2 cos(f). (5.1) 

For Az(0) <C AE, the non-rigid pendulum is driven to rotate in a direction determined by 
AE (corresponding to the ac Josephson-like effect). With AE = 0, and A > A c (of Eq. ( [4.9D , 
we had found that the pendulum also executes rotatory motion, in a direction determined 
by z(0). For Az(0) ^> AE ^ 0, we expect this type of motion to persist (corresponding 
to MQST due to nonlinearity) . In between there should be competition between the two 
effects, and a transition at some shifted critical value A = A C (AE). This physical picture 
for AE is confirmed by obtaining z(t) in terms of Weierstrassian elliptic functions that 
change their behavior at a singular value A = A C (AE). 

We show in Figs. ([)]) that the MQST phenomena, (inverse-period dip, average non-zero 
imbalance), persist in the AE ^ case, and display dependence on A and AE. Figure |9| 
shows the scaled inverse period r ac /r versus the scaled nonlinearity ratio A/A C (AE) where 
r ac is as in Eq. ((Op, with z(0) = 0.1, AE = 1.0, and 0(0) = 0,7r. The dip to zero at 
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the onset of MQST is clearly seen. The inset shows the time-averaged (z) for 0(0) = 0, it, 
vanishing at A = A C (AE). Whereas for AE = and A < A C (AE = 0), the average 
population imbalance was zero, for AE ^ we have (z) ^ in the corresponding sub- 
critical region A < A C (AE). This is analogous to a voltage across a capacitor inducing 
a charge difference and the external static magnetic field in the case of 3 He-A Note that 
there is a combined influence of A, AE, 0(0), so (z) can be larger (in magnitude) than z(0). 



In particular, for A -> 0, (z) -> -AE(^l - z 2 (0) cos 0(0) - AEz(0))/(l + AE 2 ), that for 
0(0) = is negative, as in the inset. This corresponds to an averaged pendulum rotation 
(z) ~ —AE < 0, opposite in sign to the initial z(0) > 0, but slowing to zero as the critical 
value is approached. For A > A C (AE), in the MQST regime, the averaged rotation (z) > 
is in the initial direction of z(0) > 0, with (z) approaching the initial z(0) value for large A, 
as in the AE = case of Fig. [8[ 

Shapiro Effect Analogs. Let us now consider the BJJ analog of the Shapiro resonance 
effect observed in SJJ ||23| . In addition to a time-independent trap asymmetry AE, we 
impose a sinusoidal variation so that we can write the asymmetry term as AE + AE\ cos u>ot. 
This could be done by varying the laser barrier position at fixed intensity. A similar Shapiro- 
like resonance effect could be seen, with an oscillation of the laser beam intensity, at fixed 
mid-position, so K — > /C(l + Sq coscuot). The analog of the Shapiro effect arises when the 
period from the time independent asymmetry, ~ 1/ AE, matches that from the oscillatory 
increment, ~ This matching condition is intimately connected with the phenomenon of 

Bloch oscillations and dynamic localization in crystals and trapping in two-level atoms p8| 



The dc value of the drift current, (z(t)), as a function of AE, will show up as resonant spikes. 
(For SJJ, with current drives, the Shapiro effect shows up as steps in the I-V characteristics.) 
Of course, the dc drift cannot persist indefinitely, because the phase difference between the 
condensates on the two parts of the BJJ will cease to be a well-defined quantity once the 
population in one well drops below N min . 

Figure |Tt] shows I fa oc (z(t)) obtained from time averaging the numerical solution, with 
a small ac drive and AE ^ 0. It is plotted as a function of AE/ujq for increasing values 
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of the nonlinearity ratio A. The initial conditions are z(0) ~ = 0.045,0(0) = ir/2, for 
which A c ~ 1000 (in the absence of AE and ac driving). When A is zero, sharp peaks in 

Idc occur at the usual 'Shapiro' condition values AE oc nuo,n — 1,2, As A increases, 

however, two things happen. Firstly, multiple peaks also occur at AE/u values different 
from integers. Close to the MQST regime, (A ~ A c ), there is a proliferation of peaks as 
the system moves from a regime of constant current, (z) ^ (A small), to one of constant 
population imbalance (z) ^ (A large). Secondly, the magnitude of the peaks or dc currents 
decreases. 

Finally, we note that for AE larger than the Bogoliubov quasiparticle gap A qp , and high 
enough temperatures, a dissipative quasiparticle branch might be observable. 

VI. SUMMARY 

We have investigated the Josephson dynamics in two weakly linked Bose-Einstein con- 
densates forming a Boson Josephson junction. In the resulting nonlinear two-mode model, 
we have described the temporal oscillations of the population imbalance of the conden- 
sates in terms of elliptic functions. Our predictions include non-sinusoidal generalizations of 
Josephson 'dc', 'ac' and Shapiro effects. We also predict macroscopic quantum self-trapping 
which is a self-maintained population imbalance across the junction due to atomic self- 
interaction, and 7r-oscillations, in which the phase difference across the junction oscillates 
around ir. We clarify the connection and the differences between these phenomena and oth- 
ers occuring in related systems like the superconducting Josephson junctions, the internal 
Josephson effect in 3 He-A, and Josephson oscillations between two weakly linked reservoirs 
of 3 He-B. Through a set of functional relations, we also predict the collapse of experimental 
data (corresponding to different trap geometries and total number of condensate atoms) 
onto a single universal curve. These effects constitute experimentally testable signatures of 
quantum phase coherence and the superfluid character of weakly interacting Bose-Einstein 
condensates. 
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APPENDIX A: MICROSCOPIC DERIVATION OF THE BOSON JOSEPHSON 
EQUATION FROM THE GROSS-PITAEVSKII EQUATION 



The values of the constant parameters in BJJ Eqs. ( |2.3| ) /C, E , and U depend on the 
geometry (and effective dimensionality) of the system and the total number of condensate 
atoms. We now outline their dependence in term of spatial GPE wave functions, elucidating 
the approximations underlying the BJJ equations. 



We look for the solution of the (time-dependent) GPE (|2.1|) with the following 

variational ansatz: 

tf(r,t) =M t )Mr) + M t )M r ) ( Al ) 
There are two approximations underlying this ansatz: 

1. We describe the temporal evolution of the GP wave function as the superposition of 
two wave functions (roughly) describing the condensate in each trap. The nonlinear 
interaction in GPE destroys such superposition. In effect if the condensate density in 
the tunneling region is small (as it is the case for weak links) nonlinear interaction in 
that region is negligible, and the superposition ansatz is preserved. 

2. We factorize the temporal and the spatial dependence of the GPE wave function 
describing the condensate in each trap. Later in this section, we will discuss the limit 
of validity of this approximation. 

The spatial dependence of $1,2 (r) can be constructed by the exact symmetric $+(r) and 
antisymmetric $_(r) stationary eigenstates of the GPE (see Section |TV|): 

*i(r) = (A2a) 
Mr) = (A2b) 
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ensuring that: 

J $i(r)$ 2 (r)rfr = (A3) 
and where we impose the normalization condition: 

J \<$> 1)2 {r)\ 2 dr = 1 (A4) 



Replacing Eqs. ( [Al] , [A~2] ) in the GPE ( [2.1|) , and using the orthogonality condition Eq. ( [A 3D , 
we obtain the BJJ equations, 

ih^- = (E® + Z7iJVi)^ - Ka^ (A5a) 

= {E° 2 + U 2 N 2 )ij 2 - K^ x , (A5b) 

with constant parameters : 

r h 2 

E = J ^|V^| 2 + m 2 V ext (r)dr (A6a) 

U = g J \^\ A dr (A6b) 



K 



r t2 



— (WiW 2 ) + ViV ext V2 
2m 



dr (A6c) 



We now return to our variational ansatz \I/(r, t) = ipi(t)§i(r) + ip 2 (t)& 2 (r) . The parameters 
U,AE ~ wave-function overlaps are Ay-dependent, but are independent of z(t), so the 
chemical potential difference is considered linear in z. This approximation captures the 
dominant ^-dependence of the tunneling equations coming from the scale factors ipi 2 oc 



Nx,2i but ignores shape changes in the wavefunctions for N±(t) ^ N 2 (t). We can estimate 
such corrections to the chemical potential difference A/i = fx\ — fi 2 , within the Thomas- 
Fermi approximation Hx,2 ~ N 2 ^ ~ (Ar/2) 2 / 5 (l ± z) 2 ^ 5 . Then relative corrections to the 
linear form A/i = (4/5)2 are estimated by £ = (Afi(z) — Az/5)/Afi(z) where Afi(z) = 
(1 + z) 2 / 5 - (1 - z) 2/5 . We find that £ is negligible over the z range where MQST effects 
are expected: £ ~ 0.1% for z = 0.1 and £ ~ 3% for z = 0.4. Thus Eq. QT6D with AE,A 
treated as constants, is indeed a reliable nonlinear equation describing BJJ dynamics for a 
large range of z(t) values. Similar conclusions has been reached in As a further test, 
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the GPE ( |2.1| ) has been solved numerically (in a spatial grid) in the double well geometry 
p9| , fully confirming the conclusions just outlined. 



APPENDIX B: EXACT SOLUTIONS IN TERMS OF JACOBIAN ELLIPTIC 

FUNCTIONS 



The total energy of the system is given by 

Az 2 



AEz -Vl-z 2 cos (j) = H(z{0), 0(0)) = H 0l 



(Bl) 



where H is the initial (and conserved) energy Combining Eqs. (|2.6a|) , (|2.8| ), we have 



z 2 + 



Az 2 



1 2 



+ AEz - H 



1-z 2 



(B2) 



The nonlinear GPE tunneling equations for the macroscopic amplitudes ipi^)? ^(t) are for- 
mally identical to equations governing a physically very different problem — a single electron 



in a polarizable medium, forming a polaron |32[|. Solutions have been found |32|J40| -|42|| , for 
the discrete nonlinear Schrodinger equation (DNLSE) describing the motion of the polaron 
between two sites of a dimer. Similarly we use Eq. ( |B2"D to obtain the exact solution for 
z(t) in terms of quadratures, 

dz 



At 

T 



*(0) 
z(t) 



GO a 



z 2 ) - 



2 _i_ 2zAE _ 2Hq 
'A A 



(B3) 



We consider AE = 0, and AE ^ cases separately. For symmetric double wells, AE = 0, 
the denominator of Eq. ( p3|) can be factorized, so 



At 

T 



z(0) 



dz 



*(*) J(a 2 + z 2 )(C 2 -z 2 ) 



(B4) 



where 



2 

A2 



^oA - 1) + 



; a 



A 2 



e - (h a - 1) 



C 2 (A) = 2y/A 2 + 1 - 2H A. 



(B5a) 
(B5b) 
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The solution to Eq. ( p4|) is written in terms of the 'cn' and 'dn' Jacobian elliptic functions, 
(with k, the elliptic modulus [f43f)as 



z(t) = Cca[(CA/k)(t - t ), k) for < k < 1 

= Cdn[(CA)(t - t ), l/k] for k > 1; (B6a) 

(B6b) 



1 



2 _ 1 j^J^A y 
- 2 U(A) J 2 



! + (#oA - 1) 



VA 2 + 1 - 2# A 



t = 2[AVC 2 + a 2 F(&rccos(z(0)/C), k)}' 1 , (B6c) 

where F(<j), k) = Jq d(j){l — k 2 sin 2 </>)~ 1//2 is the incomplete elliptic integral of the first kind. 

The Jacobian elliptic functions cn(w, k) and dn(u, k) are periodic in the argument u 
with period 4K(k) and 2K(k) respectively where K(k) = F(tt/2, k) is the complete elliptic 
integral of the first kind. The character of the solution changes when elliptic modulus 
k — 1. From Eq. ( P6b|) , this mathematical condition or singular parameter-dependence of 
the elliptic functions corresponds to the physical condition H = 1, A = A c of Eq. Q4.9| ), for 
the onset of MQST: fc 2 (A c ) = 1. When k 2 <C 1, cn(w, k) w cosm + fc 2 sinu(M — |sin2u) is 
almost sinusoidal. When k 2 increases, the departure from simple sinusoidal forms becomes 
drastic. For k 2 ^ 1, cn(/u, k) ~ sechw — ^-^-(tanhMsechM)(sinhM cosh-u — u) becomes non- 
periodic. When k 2 3> 1, the behavior is again periodic (but about a non-zero average): 
dn(u, l/k) w 1 - (sin 2 w)/2A; 2 . 

The time-period of oscillation of z(t) is given |E| by 

AkK(k) r . . 

r = ; for < Jfc < 1, B7a 

CA v ' 

= ?^m imk>1 . ( B7b) 

CA v 7 

In the linear limit, r — > n/y/l + A in agreement with the expression for r p in Eq. ( [4.71) . 
As A; — > 1, or A — > A c , the period becomes infinite, as in 'critical slowing down', diverging 



logarithmically, K(k) — > log(4/vl — A; 2 ). The evolution of the imbalance is given, in this 
special case, by the non-oscillatory hyperbolic secant (C = 2^A C — 1/A C ): 

z(t) = Ccn[(CA c )(t - t ), 1] = CsechCA c (t - t ) for k = 1. (B8) 
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We now turn to the AE ^ case. The general form of the integral of Eq. (p3|) is split into 
two parts, 



A/ ^o + r-d£ = 
Jm ftp) 



2 2 



where At /2 is the integral from z\ to ^(0), and Z\ is a root of the quartic f{z) = (j^j (1 — 

r 1 2 

z 2 ) — z 2 + 2z ^ E — 2 -^a . Taylor expanding f(z) around z, and with the change of variable 
y = y(z) = (f'(z 1 )/A)(z - zi)' 1 + f"(z l )/2A for which y(zi) = oo, the integral in Eq. ( |B9[) 
is cast in a standard form 

A(t-to) _ r°° rfy' , . 

2 7, y/ltf* - g 2 y' - 1 J 

which can be inverted as a Weierstrassian elliptic function y = p(A(t — to)/2; g 2 , cfe). Thus 

z(f] _ , , /'(*i)/ 4 (Bll) 

[ ) " 1 + rfA(f - t )/2; ^3) - J"(^x)/24 



In Eq. ( BIO ), the constants in the cubic h(y) = Ay 3 — g 2 y — #3 are determined from the 



coefficients a-j of f(z) = Z)f=o a «+4 2;i as 

g 2 = —04 — 4aia3 + ; #3 = —0204 + 2aia2£t3 — a 2 + — afa^, (B12) 

where 

2 , , , . q . 2# A£ 4(1 - # n 2 ) 

ai = — ! aa = ^( A (#° + 1) - AE ) ; a 3 = —jL— ; a 4 = 1 A2 ° J , (B13) 



The solution Eq. ( pil|) is equivalent to that found in polaronic j|4| and other contexts 



45| , |46[| , where AE 7^ corresponds to a difference or 'disorder' in on-site electronic or 
excitonic energies. 

For AE = we found that the elliptic modulus k 2 governed the behavior of the Jacobian 
elliptic function solutions. For AE 7^ 0, the discriminant 8 — g\ — 27g 2 . of the cubic h(y) 



(with roots 2/1,2,3) governs the behavior of the Weierstrassian elliptic functions |43|] For 5^0, 
the solutions are oscillatory about a non-zero average, (z) 7^ 0. For 5 = 0, (z) = 0, and the 
time-period diverges, corresponding to A = A C (AE), the onset of MQST. 
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The time-period of oscillation can be written in terms of complete elliptic integrals of the 
first kind K(k) as in the AE = case of Eq. (|B~7|). However, the argument and prefactors 
are different, with 

r = K(k 1 )/{y 1 - y 3 ), for 5 > 0, (B14a) 
r = K(k 2 )/^H~ 2 M 5 < 0, (B14b) 
T = oo,for5 = 0,£ 3 <0. (B14c) 

For 5 > 0, k\ = (1/2 — I/3)/ (2/1-2/3) where the roots yi of h(y) are all real, yi = — J~g2~/3 cos ([9+ 
2vr(i - l)]/3) and 9 = aiccos(yj27g$/g$). For 5 = -\5\ < 0, k 2 = 1/2 - 3j/ 2 /4(3j/| - # 2 ) where 
j/2 is the only real root, y 2 = [(g 3 + ^J-5/27) 1 ^ + (5-3 - -5 /27) 1/3 }/2. Thus the inverse 
oscillation period 1/r, for 5 7^ 0, is obtained as above in terms of A and AE, with 1/r = 
at A = A C (A£). 
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FIGURES 

FIG. 1. The asymmetric double well trap for two Bose-Einstein condensates with JV12 and 
E®2, the number of particles and the zero-point energies in the traps 1 and 2 respectively. 

FIG. 2. Population imbalance z(t) as a function of dimensionless time 2JCt (in units of h), with 
conditions A = 10, 0(0) = 0, symmetric trap. The initial population imbalance z(0) takes the 
values: a)0.1, b)0.5, c)0.59, d)0.6, e)0.65. 

FIG. 3. z(t) as a function of 2/Ct with initial conditions z(0) = 0.6, (f)(0) = tt, symmetric trap. 
A takes the values: a)0.1, b)l.l, c)l.lll, d)1.2, e)1.25, f)1.3. 

FIG. 4. The ^-potential V(4>) (in arb. units) plotted against 0/V for A = 0.2,0.4,0.6. 

FIG. 5. The z-potential W(-z) (in arb. units) plotted against z in a) and c) and the correspond- 
ing 0-evolution shown in b) and d). In a),b) 0(0) = 0, and in c),d) 0(0) = tt. The values of z(0) 
are as shown. 

FIG. 6. The inverse period (scaled in units of 2JC) 1/r plotted against z(0)/z c for a)0(O) = 0, 
b)0(O) = tt. In a), the dashed line corresponds to A = 10, for which the dip occurs at z c = 0.6 and 
the solid line to A = 100, for which z c = 0.2. The inset in a) shows the time-averaged population 
imbalance (z) as a function of z(0)/z c . In b), A = 1.111, z c = 0.6. 

FIG. 7. The phase-plane portrait of the dynamical variables z and for A/A c values as 
marked. For all trajectories, z(0) = 0.6. See text for explanation of the markings of the various 
trajectories. 

FIG. 8. Universal curve for data-collapse with (tt / '(ICC A 2 t) (in units of Ti) versus A; 2 (A) as in 
Eq. (p6b| ). 
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FIG. 9. Scaled inverse period r ac /r plotted against A for fixed asymmetric trap parameter 
AE = 1, z(0) = 0.1, and 0(0) = 0, it initial values, and l/r ac as defined in Eq. ( [4.7|) . The vertical 
scale on the left (right) corresponds to 0(0) = 0(7r). The insets show time-averaged {z) against A, 
for (a), 0(0) = and (b), 0(0) = vr. 

FIG. 10. Analog of Shapiro effect: dc current I^c = (i) versus trap asymmetry parameter 
scaled in applied frequency, AE/ojq. Here z(0) = 0.045,0(0) = w/2, AEi/ujqTi = 3.5, and dashed 
(thick solid) lines are for A = 0(1000). 
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